Shrinking of a condensed fermionic cloud in a trap approaching the BEC limit 
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We determine the zero-temperature density profile of a cloud of fermionic atoms in a trap subject 
to a mutual attractive interaction, as the strength of the interaction is progressively increased. We 
find a significant decrease of the size of the atomic cloud as it evolves from the weak-coupling (BCS) 
regime of overlapping Cooper pairs to the strong-coupling (Bose-Einstein) regime of non-overlapping 
bound-fermion pairs. Most significantly, we find a pronounced increase of the value of the density 
at the center of the trap (even by an order of magnitude) when evolving between the two regimes. 
Our results are based on a generalized Thomas-Fermi approximation for the superfluid state, that 
covers continuously all coupling regimes. 
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The impressive experimental realization of Bosc- 
Einstcin condensation (BEC) of ultracold bosonic atoms 
in a trap over the last few years 1 has prompted the ex- 
pectation that also the properties of ultracold fermionic 
atoms in a trap could be revealed in a regime of quan- 
tum degeneracy 2 . In addition, the possibility that the 
study of Fermi gases may lead to a superfluid phase anal- 
ogous to the BCS state has added excitement to the field. 
In particular, experimental proposals are under investi- 
gation for realizing Cooper pairs with ultracold dilute 
fermionic atoms in a trap, either via the enhanced inter- 
action between fermionic atoms with a tunable scattering 
(Feshbach) resonance 3 or via sympathetic cooling with a 
mixture of fermions and bosons 4 . 

Given these rapid experimental advances, it seems 
timely to explore theoretically the even more in- 
triguing question of the crossover between the weak- 
coupling (BCS) regime of overlapping Cooper pairs 
and the strong-coupling (Bose-Einstein) regime of non- 
overlapping bound-fermion pairs for a gas of ultracold 
Fermi atoms in a trap, thus dealing simultaneously with 
two basic quantum phenomena that share the same spon- 
taneous broken-symmetry (superfluid) behavior. In par- 
ticular, novel and interesting physics may result in the 
intermediate-coupling regime where neither one of the 
two limits is fully realized. 

Detection of density profiles of ultracold trapped atoms 
is routinely used to explore their degeneracy properties. 
This feature could also be explored to detect the con- 
tinuous evolution from Cooper pairs to bound-fermion 
pairs as the strength of the attractive interaction between 
fermions is increased, with the expectation that the den- 
sity profile shrinks upon approaching the strong-coupling 
regime owing to the diminished Fermi pressure. In this 
respect, partial results have already been obtained as a 
function of temperature in the broken-symmetry phase 
for a single coupling toward the weak-coupling side of the 
BCS-BEC crossover 5 , as well as a function of coupling in 



the normal phase above the critical temperature 6 . 

In this paper, we calculate the density profile of a cloud 
of fermionic atoms in a trap subject to a mutual attrac- 
tive interaction, as its strength is progressively increased 
from weak to strong coupling. To this end, we extend 
the coupled mean-field equations determining the gap 
function and the chemical potential for the BCS-BEC 
crossover problem of a homogeneous system to the case 
of a spatially varying gap function in the presence of a 
harmonic trap, and study in a systematic way the evolu- 
tion of the density profile at zero temperature where the 
mean-field equations are expected to be reliable. Our 
approach generalizes the fermionic Thomas-Fermi (TF) 
approximation to the superfluid state at zero tempera- 
ture for any coupling, and recovers the TF approxima- 
tion for the composite bosons in the strong-coupling 
limit. We find a significant decrease of the size of the 
atomic cloud as the strong-coupling limit is approached; 
even more significantly, we find a pronounced increase of 
the value of the density at the center of the trap when 
evolving between the two regimes. This increase can thus 
be regarded as a signature of the achieved BCS-BEC 
crossover. Numerical calculations for the density profile, 
gap function, and chemical potential, as well as analytic 
results for the radius of the cloud and related quantities 
are presented. 

A central issue of the BCS-BEC crossover problem 8 
is to obtain for all coupling regimes a single theory ex- 
pressed in terms of fermionic variables, which recovers 
known results for weakly-interacting fermions (on the 
weak-coupling side) and for weakly interacting bosons 
(on the strong-coupling side), thus providing an interpo- 
lation scheme over the intermediate coupling regime. The 
treatment of the intermediate-coupling regime is neces- 
sarily an approximate one, because in this regime a small 
parameter to control the approximations is lacking. For 
these reasons, the interpolation via a single fermionic the- 
ory represents at present the only manageable approach 
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to the BCS-BEC crossover. [Effects going beyond the 
BCS theory in weak coupling have been included for the 
homogeneous case in Ref. 9 In the present treatment of 
the inhomogeneous case, we shall limit to consider the 
crossover from BCS theory to the BEC limit, thus leav- 
ing out these additional effects.] 

The present calculations are based on a local-density 
approximation, whereby the fermionic chemical poten- 
tial entering the BCS gap and density equations is re- 
placed by a local chemical potential that adapts to the 
spatial position in the harmonic trap. Wc thus set 
yu(r) = /i — V(r), where /x is the equilibrium chemical 
potential and V(r) — mw 2 r 2 /2 is the harmonic potential 
for the Fermi atoms at position r in the trap (m being the 
fermionic mass and w the characteristic trap frequency). 
Correspondingly, the gap function takes the local value 
A(r). With the replacement of ti by (i(r) and of A by 
A(r), the s-wave BCS gap equation at zero temperaturen 
becomes r dependent and reads in dimensionless form: 
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(1) 

We have here introduced the dimensionless variables 
f = \t\/R f , A(f) - A(\r\)/E F , and /2(f) = (i(\r\)/E F = 
Ji — f 2 for a spherical trap, where E F = (SN) 1 / 3 ^ (h = 
1 throughout) is the Fermi energy for non-interacting 
fermions in the trap 7 , R F is identified by the condition 
moj 2 R F /2 = E F , and fi = (i/E F . We have further de- 
fined the Fermi wave vector k F by E F = hp /(2m). 

The above equation has been obtained by considering a 
contact potential for the attractive fermionic interaction 
active between two equally populated fermionic species 
with different internal ( "spin" ) states, and then regular- 
izing the ensuing ultraviolet divergence of the gap equa- 
tion in terms of the scattering length a F of the associated 
fermionic two-body problem 10 ' 11 . This regularization, in 
fact, considerably simplifies dealing with the BCS-BEC 
crossover and parametrizes the interaction in terms of the 
measurable scattering length a F . 

Dealing with the BCS-BEC crossover requires us to 
supplement the gap equation by the density equation 
N = Jdrn(r) = Jdrh(r), yielding 
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Here, f c is the radius of the atomic cloud defined by 
the condition n(f) = for r > f c (which, in turn, re- 
quires A(f) = for f > r c as a necessary condition). 
It can be readily verified that f c is determined by the 
equation (i(f c ) — for a F < 0, yielding f c — \f]l, 
and by the equation /2b (f c ) = for a F > 0, yielding 
r c = (k F a F )~ 2 . In the above expression, we have 



introduced the local bosonic chemical potential /Us(r) = 
fiB — Vb(y), where (1b = 2/i + eq is the equilibrium 
bosonic chemical potential, Vg(r) = 2V(r) = m B u 2 r 2 /2 
is the harmonic potential for the bosonic molecules of 
mass niB — 2m, and eq = (ma 2 F )~ l is the molecular 
binding energy. Contrary to Eq. (1), Eq. (2) couples 
different positions in the trap, since obtaining a uniform 
chemical potential (i requires knowledge of A(r) and /i(r) 
over the whole trap region. 

The weak- and strong-coupling limits of Eqs. (1) 
and (2) are essentially determined by the behavior of 
the chemical potential (i, which about coincides with the 
Fermi energy E F in the weak-coupling limit and with 
(minus half) the binding energy e in the strong-coupling 
limit. These values represent, in fact, the energy required 
to extract one fermion from the system in the two limits, 
namely, from the Fermi sea and from a bound molecule, 
respectively. 

In analogy to the homogeneous case 12 , analytic re- 
sults can be obtained from the coupled equations (1) 
and (2) both in the weak- and strong-coupling limits, 
defined, in the order, by the conditions (i(f) > and 
A(f) < jl(f), and (i(f) < and A(f) < |/l(f)| (for 
< f < r c ). In weak coupling, one obtains A(f) = 

(8/2(f)/e 2 ) exp j— n/ (2k F \a F \y / (l(r))^ with the typical 
essential singularity behavior as fcp|ai?| — > 0, and the 
associated density profile n(f) = (8N/n 2 ) (l — f 2 ) 3 ^ 2 
which is the standard result within the TF approximation 
for non-interacting fermions. For the radius of the atomic 
cloud one gets correspondingly r c — R F . In strong 
coupling, one obtains instead A(f) 2 = 2(eo/ E F )fj,B(r) 
with the associated bosonic density ns(r) = n(r)/2 = 
(ms/47ras)/iB(r) in terms of the original variables, 
where as = 2a F identifies the bosonic scattering length. 
The TF result for bosons 7 is thus properly recovered by 
our calculation in the strong-coupling limit, with a resid- 
ual bosonic interaction Uq = Anas /tub that matches pre- 
viuos estimates for the homogeneous case. [A more so- 
phisticated theory for the homogeneous case 11 has shown 
that the correct value for the bosonic scattering length 
as is smaller by about a factor of 3 than the value 
2a F . Owing to the very slow dependence of the size of 
the Thomas- Fermi density profile on as in strong cou- 
pling (see below), this difference makes practically no 
change in this quantity.] Using eventually the expres- 
sion of (Mb for bosons in a trap 7 , one gets the value 

= \J Mb/2 = 0.69(k F a F ) 1 ^ 5 showing a nontrivial de- 
pendence of f c on the coupling parameter in the strong- 
coupling limit. 

Recovering from Eqs. (1) and (2) the TF results both 
for non-interacting fermions (in the weak-coupling-limit) 
and for weakly-interacting bosons (in the strong-coupling 
limit) should not be surprising, since these equations can 
also be obtained by generalizing the TF procedure for 
fermions to the superfluid state. To this end, one intro- 
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duces an energy functional E[n(r), A(r)] that depends 
on the local variables n(r) and A(r), such that minimiz- 
ing E[n(r), A(r)] with respect to variations of A(r) and 
n(r) (subject to the constraint TV — Jrfrn(r)) leads to 
Eqs. (1) and (2), respectively. Within the present ap- 
proach, this functional is given by the expression 
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where £k = fc 2 / (2m) and A is here an implicit function of 
n(r) and A(r) via the equation 
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In particular, in the strong-coupling limit (whereby 
n(r) oc A(r) 2 ) the functional (3) reduces to 
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from which the bosonic TF expression for ne(r) can be 
obtained by minimizing with respect to n(r) and taking 
the constraint into account. [The absence of the kinetic 
energy term in the bosonic TF approximation results in 
the present approach by the neglecting of the center-of- 
mass motion of the pairs in the expression (3) for the 
energy.] This results implements at the level of the TF 
approximation the density functional theory for super- 
conductors described in Ref. 13. 

The rescaled equations (1) and (2) depend only on the 
dimcnsionlcss coupling parameter (kpap)^ 1 and not on 
the total number of particles TV. This is due to the fact 
that, as expected, the kinetic energy is not properly taken 
into account in the present local-density treatment that 
reduces to the TF approximation for fermions in weak 
coupling and for bosons in strong coupling. What is miss- 
ing in our equations is the presence of the length scale of 
the harmonic trap, which for a bosonic molecule equals 
a osc = (tjb^) -1 / 2 and to which there would be asso- 
ciated the TV-dependent dimensionless ratio Rp/a^ ac = 
2.402 TV 1 / 6 . The physical lower bound r c = af sc , how- 
ever, can never be reached in practice, since it would cor- 
respond to the condition (kpap)^ 1 = 12.5 TV 5 / 6 , yielding 
the totally unrealistic value a F ~ 10 _3 ao (a being the 
Bohr radius) in the case of a fermionic atom like 40 K for 
TV <~ 10 6 and a typical trap value uj ~ 10 3 s _1 . 

By a similar token, it can be shown that the local- 
density approximation for the trapped Fermi gas holds 
provided kpRp 3> 1 (such that the energy quantiza- 
tion in the trap is irrelevant with respect to Ep) and 
A(r = 0)/u » 1 (such that the pairs are well contained 
within the trap) 14 . When TV ~ 10 6 , k F R F ~ 10 2 while 



A(r = 0)/w ~ 10 already for the weak-coupling value 
(kpap)- 1 = -2. 

Figure 1 shows the density profile n(r) along the ra- 
dial coordinate (in units of R F ) for several values of the 
coupling parameter {kpap)^ 1 . The TF profiles for non- 
interacting fermions as well as for interacting point-like 
bosons (with (fcj?ap) _1 = 5) are also shown for com- 
parison. Note from the figure that the non-interacting 
fermion density profile is recovered in weak coupling 
and the interacting point-like boson density profile is re- 
covered in strong coupling, respectively. Note also the 
shrinking of the size of the atomic cloud as it evolves 
from the weak-coupling to the strong-coupling regime, 
the critical radius f c being reduced by about a factor of 
2 across the crossover region. Correspondingly, the den- 
sity at the center of the trap increases significantly (by 
about a factor of 6 for the coupling values shown in the 
figure) . 
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FIG. 1. Density profile (divided by TV) vs the radial coor- 
dinate for several values of the coupling parameter (fcffli?) -1 . 

Figure 2 shows the radius f c of the atomic cloud (full 
line) vs the coupling parameter, as determined by enter- 
ing the numerical values of the chemical potential (ob- 
tained from the solution of Eqs. (1) and (2)) into the ana- 
lytic expressions of r c reported previously. Three distinct 
regions of the coupling parameter are identified from 
this plot: (i) A weak-coupling region for (kpap)^ 1 < 
— 1 where r c is almost coupling-independent and equal 
to 1; (ii) An intermediate-coupling region for —1 < 
{kpap)^ 1 < 1 where r c rapidly decreases for increasing 
coupling, being reduced by about 35% when (kpap)^ 1 = 
1; (iii) A strong-coupling region for 1 < (kpap)^ 1 where 
f c coincides with the analytic expression of the bosonic 
TF theory (dashed line). For completeness, we also show 
in the inset the value of the density n(f — 0) at the center 
of the trap vs the coupling parameter (full line). Follow- 
ing a flat behavior in the weak-coupling regime, a quite 
rapid increase of h(f = 0) occurs in the intermediate- 
coupling regime, approaching eventually the power-law 
dependence n(f = 0) cx (kpap)~ 3 / 5 obtained analyti- 
cally in the strong-coupling regime (dashed line). 
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The chemical potential is finally shown in Fig. 3 (full 
line) vs the coupling parameter. The behavior of the 
chemical potential obtained from the expression /j, — 
{fJ-B — £o)/2 in strong coupling (with the value of \ib for 
bosons in a trap 7 ) is also shown for comparison (dashed 
line). 
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FIG. 2. Radius of the atomic cloud (full line) and its 
strong-coupling approximation (dashed line) vs the coupling 
parameter (AfOf) -1 - Inset: corresponding values of the di- 
mensionless density at the center of the trap. 
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FIG. 3. Chemical potential in units of Ep when [i > and 
of eo /2 when fj, < (full line) and its strong-coupling approx- 
imation (dashed line) vs the coupling parameter (^fchf) -1 • 
Inset: gap parameter at the center of the trap (full line) 
and its strong-coupling approximation (dashed line) vs the 
coupling parameter (the value for 3D homogeneous case 
(dashed-dotted line) is also shown for comparison). 

We are indebted to M. Inguscio and G. Modugno for 
discussions. 



In the inset, the gap parameter A(f = 0) at the center 
of the trap (full line) is shown vs the coupling param- 
eter and compared with its strong-coupling approxima- 
tion (dashed line) (the value for the three-dimensional 
homogeneous case 12 (dashed-dotted line) is also shown). 
Note the marked difference between the trapped and ho- 
mogeneous cases, with the value for the trapped case 
increasing much faster than for the homogeneous case 
as the coupling increases. Specifically, in strong cou- 
pling the analytical arguments discussed above yield 
A(f) = {2y / 2/k F a F ) y / 0A76(k F a F ) 2 / 5 - P, producing 
the value A(0) = 1.951(k F a F )- 4 / 5 at the center of the 
trap; for the homogeneous case, on the other hand, 
A/E F = 1.303(fc F a F )- 1 / 2 . 

As the data of Figs. 1-3 for the various physical quan- 
tities n(r), A(r), r c , and n show, the BCS-BEC crossover 
takes place in practice over the rather narrow range 
— 1 < (kpap)^ 1 < 1 of the coupling parameter. This re- 
mark suggests that it would be especially interesting to 
explore experimentally this intermediate-coupling region, 
where deviations from the purely fermionic and bosonic 
behaviors occur. 

In conclusion, we have shown that the density profile 
for superfluid fermionic atoms in a trap, exhibits an in- 
teresting shrinking of the size of the atomic cloud as the 
strength of the attraction increases. This behavior should 
permit one to decide whether the superfluid behavior is 
either of the BCS type with strongly overlapping Cooper 
pairs or of the BEC type with non-overlapping bound- 
fcrmion pairs. 
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